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Abstract. We show that vector theories on cosmological scales are excellent candidates for dark 
energy. We consider two different examples, both are theories with no dimensional parameters 
nor potential terms, with natural initial conditions in the early universe and the same number of 
free parameters as ACDM. The first one exhibits scaling behaviour during radiation and a strong 
phantom phase today, ending in a "big-freeze" singularity. This model provides the best fit to date 
for the SNIa Gold dataset. The second theory we consider is standard electromagnetism. We show 
that a temporal electromagnetic field on cosmological scales generates an effective cosmological 
constant and that primordial electromagnetic quantum fluctuations produced during electroweak 
scale inflation could naturally explain, not only the presence of this field, but also the measured 
value of the dark energy density. The theory is compatible with all the local gravity tests, and is 
free from classical or quantum instabilities. Thus, not only the true nature of dark energy could be 
established without resorting to new physics, but also the value of the cosmological constant would 
find a natural explanation in the context of standard inflationary cosmology. 
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INTRODUCTION 

The fact that today matter and dark energy have comparable contributions to the energy 
density, ~ Pm ~ (2 x 10^-^ eV)^ in natural units, poses one of the most important 
problems for models of dark energy. Thus, if dark energy is a cosmological constant, its 
energy density would remain constant throughout the history of the universe, whereas 
those of the rest of components (matter and radiation) grow as we go back in time. Then 
the question arises as to whether it is a coincidence (or not) that they have comparable 
values today when they have differed by many orders of magnitude in the past. Notice 
also that if A is a fundamental constant of nature, its scale (around 10^^ eV) is more 
than 30 orders of magnitude smaller than the natural scale of gravitation, G — Mp 
with Mp ~ 10^^ GeV. On the other hand, alternative models in which dark energy is a 
dynamical component rather than a cosmological constant also require the introduction 
of unnatural scales in their Lagrangians or initial conditions in order to account for the 
present phase of accelerated expansion. Such models are usually based on new physics, 
either in the form of new cosmological fields or modifications of Einstein's gravity 
ill, 0, HL 3, El] and they are generically plagued by additional problems such as classical 
or quantum instabilities, or inconsistencies with local gravity constraints. 

Therefore, we would like to find a description for dark energy without dimensional 
scales (apart from Newton's constant G), with the same number of free parameters as 
ACDM, with natural initial conditions, with good fits to observations and no consistency 



problems. In this work we will show that vector theories can do the job. With that 
purpose, we present two examples of such theories which have been recently proposed 
laH] (for other vector models see references in those works). 



SCALING VECTOR DARK ENERGY 

The action in this case reads fd]: 

S = J d'xV^(^-j^-^-F,,F^^-^-iV,A^)^+R,,A^A'^ (1) 

Notice that the theory contains no free parameters, the only dimensional scale being the 
Newton's constant. The numerical factor in front of the vector kinetic terms can be fixed 
by the field normalization. Also notice that the "mass" termT^^y^'^^^ can be written as 
a combination of derivative terms as V^A'^VyA^ — V^A^VyA^ and therefore the theory 
contains no potential terms. This action resembles that of Maxwell electromagnetism in 
the Feynman gauge with a mass term. 

The classical equations of motion derived from the action in ([T]) are the Einstein's and 
vector field equations: 

Rnv-^Rguv = S7rG{T^y + T^y) (2) 
UA^+R^yA" = (3) 

where T^y is the conserved energy-momentum tensor for matter and radiation and T^y 
is the energy-momentum tensor coming from the vector field. For the simplest isotropic 
and homogeneous flat cosmologies, we assume that the spatial components of the vector 
field vanish, so that = {Ao{t), 0,0,0) and that the space-time geometry will be given 
by: 

ds^ = dt^ - a^{t) dijdx'dx^ , (4) 

For this metric ([3]) reads: 

Aq + 3HAq - 3 [2H^ +H]Aq = (5) 

Assuming that the universe has gone through radiation and matter phases in which the 
contribution from dark energy was negligible, we can easily solve this equation in those 
periods. In that case, the above equation has a growing and a decaying solution: 

Ao(0=A+r«++Aor«- (6) 

with Aq constants of integration and a± = — (l±l)/4in the radiation era, and a± = 

(—3 ± a/33)/6 in the matter era. On the other hand, the (00) component of Einstein's 
equations reads: 
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where the vector energy density is given by 
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Figure 1: (Left) Evolution of energy densities for the best fit model. Dashed (red) for radiation, dotted 
(green) for matter and solid (blue) for vector dark energy. We show also for comparison the cosmological 
constant density in dashed-dotted line. (Right) Evolution of dark energy equation of state for the best fit 
model. The lower panel shows the 1 a confidence interval. 



Using the growing mode solution from we obtain = pAoa^ with K = —4 in 
the radiation era and k = {\/33 — 9)/2 ~ —1.63 in the matter era. Thus, the energy 
density of the vector field starts scaling as radiation at early times, so that Pa/Pr = 
const. However, when the universe enters its matter era, starts growing relative to 
pM eventually overcoming it at some point, in which the dark energy vector field would 
become the dominant component (see Fig. 1). Notice that since Aq is essentially constant 
during radiation era, solutions do not depend on the precise initial time at which we 
specify it. Thus, once the present value of the Hubble parameter Hq and the constant 
Ao during radiation (which fixes the total matter density Q.m) are specified, the model 
is completely determined, i.e. this model contains the same number of parameters as 
ACDM, which is the minimum number of parameters of a cosmological model with 
dark energy. As seen from Fig. 1 the evolution of the universe ends at a finite time tend 
where a Uend with a^nd finite, Ao(w) = Mp/ (4^/^), poE °° and poE ^ This 
corresponds to a Type III (big -freeze) singularity according to the classification in [|8|]. 

We can also calculate the effective equation of state for dark energy as: 

^PA^ -3{lH^ + lH)Al+HAoAo-lAl 
"^^^ Pa lH^Al + 3HAoAo-^Al 

Again, using the approximate solutions in we obtain: wde = 1/3 in the radiation era 
and Wde — —0.457 in the matter era. As shown in Fig. 1, the equation of state can cross 
the so called phantom divide, so that we can have wi)e{z = 0) < —1. 

In order to confront the predictions of the model with observations of high-redshift 
supernovae type la, we have carried out a statistical analysis for two supernovae 
datasets, namely, the Gold set [|9|], containing 157 points with z < 1.7, and the more 
recent SNLS data set [|lO|] . comprising 115 supernovae but with lower redshifts (z< 1). 
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Table 1: Best fit parameters with la intervals for the vector model (VCDM) and the cosmological 
constant model (ACDM) for the Gold (157 SNe) and SNLS (115 SNe) data sets, wq denotes the present 
equation of state of dark energy. Aq is the constant value of the vector field component during radiation. 
Zt is the deceleration-aceleration transition redshift. fo is the age of the universe in units of the present 
Hubble time. tg„ii is the duration of the universe in the same units. 

In Table 1 we show the results for the best fit together with its corresponding lo 
intervals for the two data sets. We also show for comparison the results for a standard 
ACDM model. We see that the vector model (VCDM) fits the data considerably better 
than ACDM (in more than 2a) in the Gold set, whereas the situation is reversed in the 
SNLS set. This is just a reflection of the well-known 2a tension [fllll between the two 
data sets. Compared with ACDM, we see that VCDM favors a younger universe (in H^^ 
units) with larger matter density. In addition, the deceleration- acceleration transition 
takes place at a lower redshift in the VCDM case. The present value of the equation of 
state with wq = —3.53^^'^^ which clearly excludes the cosmological constant value —1. 



Future surveys 111211 are expected to be able to measure wq at the few percent level and 
therefore could discriminate between the two models. 

We have also compared with other parametrizations for the dark energy equation of 



state [|13n . Since our one-parameter fit has a reduced chi-squared: x /d-o.f = 1.108, 



VCDM provides the best fit to date for the Gold data set. 

We see that unlike the cosmological constant case, throughout radiation era p^E / pR ~ 
10^^ in our case. Moreover the scale of the vector field Aq = 3.71 x 10^^M/> in that era is 
relatively close to the Planck scale and could arise naturally in the early universe without 
the need of introducing extremely small parameters (for instance in an inflationary 
epoch), thus avoiding the coincidence problem. 

In order to study the model stability we have considered the evolution of met- 
ric and vector field perturbations. Thus, we obtain the dispersion relation and the 
propagation speed of scalar, vector and tensor modes. For all of them we obtain 
V = (1 — \6tzGA^^^^^ which is real throughout the universe evolution, since the value 
Aq = {l6nG)^^ exactly corresponds to that at the final singularity. Therefore the model 
does not exhibit exponential instabilities. As shown in Hill , the fact that the propagation 
speed is faster than c does not necessarily implies inconsistencies with causality. We 
have also considered the evolution of scalar perturbations in the vector field generated 
by scalar metric perturbations during matter and radiation eras, and, again, we do not 
find exponentially growing modes. 

If we are interested in extending the applicability range of the model down to so- 
lar system scales then we should study the corresponding post-Newtonian parameters 
(PPN). We can see that for the model in (flOl) . the static PPN parameters agree with those 
of General Relativity [1150, i.e. 7= /3 = 1. For the parameters associated to preferred 
frame effects we get: OCi = and a2 = SttAq/Mp where Aq is the norm of the vector 
field at the solar system scale. Current limits a2 ^ 10^^ (or a2 ^ 10^^ for static vector 
fields during solar system formation) then impose a bound Aq< 10^^(10^^)Mp. In or- 
der to determine whether such bounds conflict with the model predictions or not, we 
should know the predicted value of the field at solar system scales, which in principle 
does not need to agree with the cosmological value. Indeed, Aq will be determined by 
the mechanism that generated this field in the early universe characterized by its pri- 
mordial spectrum of perturbations, and the subsequent evolution in the formation of the 
galaxy and solar system. Another potential difficulty arising generically in vector-tensor 
models is the presence of negative energy modes for perturbations on sub-Hubble scales. 
They are known to lead to instabilities at the quantum level, but not necessarily at the 
classical level as we have shown previously. Work is in progress in order to determine 
which are the necessary conditions to avoid the presence of such states in this model. 



IS THE NATURE OF DARK ENERGY ELECTROMAGNETIC? 

In the previous section, we have seen that a generic vector field whose action resembles 
that of electromagnetism with a mass term could be a good candidate for dark energy, 
but what about standard electromagnetism?. In a very recent work it has been 
shown using the covariant (Gupta-Bleuler) formalism that it is indeed possible to explain 
cosmic acceleration from the standard Maxwell's theory. 



We start by writing the standard electromagnetic action including a gauge-fixing term 
in the presence of gravity: 



S = / d^xy 



167tG 4 2 



(10) 



The gauge-fixing term is required in order to define a consistent quantum theory for 
the electromagnetic field UM, and we will see that it plays a fundamental role on 
large scales. Still this action preserves a residual gauge symmetry + (9^0 with 

n(j) —0. Electromagnetic equations derived from this action can be written as: 

VyF^" + AV^VvA^ = (11) 

Notice that since we will be using the covariant Gupta-Bleuler formalism, we do not a 
priori impose the Lorentz condition. 

We shall first focus on the simplest case of a homogeneous electromagnetic field 
^/x — (^o(05^(0) ill ^ fl^t Robertson-Walker background. In this space-time, equations 
(fill) read: 

Ao + 3HAo + 3HAo = 

A + HA = (12) 

We can solve (fT2l) during the radiation and matter dominated epochs when the Hubble 
parameter is given hy H = p/t with p = 1/2 for radiation and p = 2/3 for matter. In 
such a case the solutions for (fT2l) are: 

Ao(0 = Ap+A^t-^P (13) 
A{t) = A+t^-P+A- (14) 

where Aq and are constants of integration. Hence, the growing mode of the temporal 
component does not depend on the epoch being always proportional to the cosmic time 
t, whereas the growing mode of the spatial component evolves as t^^^ during radiation 
and as t^^^ during matter, i.e. at late times the temporal component will dominate over 
the spatial ones. 

The energy densities of the temporal and spatial components read: 

^HWQ + 3HAoAo + ^Ai^ (15) 

(i)2 (16) 

Notice that we need A > in order to have positive energy density for Aq. In fact, 
it is possible to show that imposing canonical normalization for the corresponding 
creation and annihilation operators we get A = 1/3 [jvl]. Besides, when inserting the 
growing modes of the fields into these expressions we obtain that p^^ = p^^, p^ = P^<^ 
and V^A^ = const. Thus, the field behaves as a cosmological constant throughout the 
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evolution of the universe since its temporal component gives rise to a constant energy 
density whereas the energy density corresponding to A always decays as radiation. 
Moreover, this fact prevents the generation of a non-negligible anisotropy which could 
spoil the highly isotropic CMB radiation. Finally, when the universe is dominated by 
the electromagnetic field, both the Hubble parameter and Aq become constant (one can 
straightforwardly check that this is a solution of the complete system of equations) so 
the energy density is also constant and the electromagnetic field behaves once again as 
a cosmological constant leading therefore to a future de Sitter universe. As the observed 
fraction of energy density associated to a cosmological constant today is — 0.7, we 
obtain that the field value today must be Ao(?o) — 0.3 Mp. 

The effects of the high electric conductivity o can be introduced using the magneto- 
hydrodynamical approximation and including the current term 7, = o{doAi — (9/Ao) on 
the r.h.s. of Maxwell's equations. Notice that because of the universe electric neutrality, 
conductivity does not affect the evolution of Ao(?)- The infinite conductivity limit simply 
eliminates the growing mode of A{t) in (fT4)) . 

We still need to understand which are the appropriate initial conditions leading to the 
present value of Aq. In order to avoid the cosmic coincidence problem, such initial con- 
ditions should have been set in a natural way in the early universe. In a very interesting 
work flv], it was suggested that the present value of the dark energy density could be 
related to physics at the electroweak scale since pA ~ (M^y^/Mp)^, where Mew ~ lO'^ 
GeV. This relation offers a hint on the possible mechanism generating the initial am- 
plitude of the electromagnetic fluctuations. Indeed, we see that if such amplitude is set 
by the size of the Hubble horizon at the electroweak era, i.e. Ao{tEw)^ ~ ^ew^ ^^^^ 
the correct scale for the dark energy density is obtained. Thus, using the Friedmann 
equation, we find H^^, ~ M'^^/Mp, but according to (fTSi) . pA,, ~ H^A^ ~ const., so that 
Paq ~ H^y[/ ~ (A/|jy/A//>)^ as commented before. 

A possible implementation of this mechanism can take place during inflation. Notice 
that the typical scale of the dispersion of quantum field fluctuations on super-Hubble 
scales generated in an inflationary period is precisely set by the almost constant Hubble 
parameter during such period Hi, i.e. (Aq) ~ Hj |18|] . The correct dark energy density 
can then be naturally obtained if initial conditions for the electromagnetic fluctuations 
are set during an inflationary epoch at the scale M/ ~ Mew ■ 

Despite the fact that the background evolution in the present case is the same as 
in ACDM, the evolution of metric perturbations could be different, thus offering an 
observational way of discriminating between the two models. In fact, the evolution of 
the scalar perturbation with respect to the ACDM model gives rise to a possible 
discriminating contribution to the late-time integrated Sachs-Wolfe effect [|l9l] . The 
propagation speeds of scalar, vector and tensor perturbations are found to be real and 
equal to the speed of light, so that the theory is classically stable. We have also checked 
that the theory does not contain ghosts and it is therefore stable at the quantum level. 
On the other hand, using the explicit expressions in [15] for the vector-tensor theory of 
gravity corresponding to the action in (flOl) . it is possible to see that all the parametrized 
post-Newtonian (PPN) parameters agree with those of General Relativity, i.e. the theory 
is compatible with all the local gravity constraints for any value of the homogeneous 
background vector field [ 20.1 . 



The presence of large scale electric fields generated by inhomogeneities in the Aq field 
opens also the possibility for the generation of large scale currents which in turn could 
contribute to the presence of magnetic fields with large coherence scales. This could shed 
light on the problem of explaining the origin of cosmological magnetic fields. Work is 
in progress in this direction. 

CONCLUSIONS 

We have shown that vector theories offer a simple and accurate description of dark en- 
ergy in which the coincidence problem could be easily avoided. In our first example, the 
scaling behaviour during radiation and the natural initial conditions for the vector field 
offer a neat way around the problem. Moreover, in our second example, the presence 
of a cosmological electromagnetic field generated during inflation provides a natural 
explanation for the cosmic acceleration. This result not only offers a solution to the 
problem of establishing the true nature of dark energy, but also explains the value of 
the cosmological constant without resorting to new physics. In this scenario the fact that 
matter and dark energy densities coincide today is just a consequence of inflation taking 
place at the electroweak scale. Present and forthcoming astrophysical and cosmologi- 
cal observations will be able to discriminate these proposals from the standard ACDM 
cosmology. 
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